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In the Hermitian regime, the Berry phase is always a real number. It may be imaginary for a 
non-Hermitian system, which leads to amplitude amplification or attenuation of an evolved quantum 
state. We study the dynamics of the non-Hermitian Rice-Mele model driven by a time-dependent 
external field. The exact results show that it can have full real spectrum for any value of the field. 

Several rigorous results are presented for the Berry phase with respect to the varying field. We 
demonstrate that the Berry phase is the same complex constant for any initial state in a single sub¬ 
band. Numerical simulation indicates that the amplitude control of a state can be accomplished by 
a quasi-adiabatic process within a short time. 
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I. INTRODUCTION 

In the Hermitian regime, the Berry phase n is al¬ 
ways the real number. It may be complex for a non- 
Hermitian system, the concept of which was first intro¬ 
duced by Garrison and Wright 0 in a dissipative sys¬ 
tem. And in these years the geometric phase for quan¬ 
tum systems governed by non-Hermitian Hamiltonians 
and complex-valued geometric phase effects have gained 
considerable attention and have been studied by vari¬ 
ous authors (Ml- The existence of imagine part of the 
phase could lead to amplitude amplification or attenu¬ 
ation of an evolved quantum state. We investigate the 
dynamics of the non-Hermitian Rice-Mele model driven 
by a time-dependent external field. The exact results 
show that it can have full real spectrum for any value 
of the field and several rigorous results are presented for 
the Berry phase with respect to the varying field. We 
find that the Berry phase is the same complex constant 
for any initial state in a single sub-band. And via nu¬ 
merical simulations, the amplitude control of a state can 
be achieved by a quasi-adiabatic process within a short 
time. 

In a Hermitian quantum system, the geometric phase 
acquired during an adiabatic evolution is always real and 
can bring nothing to an evolved quantum state if only 
one eigenstate is involved. It has been shown that a 
non-Hermitian system can have real spectrum [l8l| and 
possess pec uliar phenomena. These include fast prop¬ 
agation |l9| . infinite reflection coefficient [20l - l^ . unidi¬ 
rectional transmission [2^ , transmission phase lapse [2^ , 
maximum multi-particle entanglement associated with 
the phase transition [2^ as well as the complex Berry 
phase. The imaginary part of the Berry phase is signif¬ 
icant for a propagating particle since it may be utilized 
to directly amplify or attenuate the particle probability. 
Very recently, the spectral and dynamical properties of 
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a quantum particle constrained on a ring threaded by a 
time-varying magnetic flux in the presence of a complex 
(non-Hermitian) potential are investigated [1^. It has 
been demonstrated that several striking effects are ob¬ 
served in the non-Hermitian case in comparison with the 
Hermitian one. 


In a previous work 27|, the dynamical behavior has 
been investigated for a non-Hermitian Rice-Mele model 
in the absence of an external magnetic field. It has 
been shown that within the unbroken V'T -symmetric re¬ 
gion, the translational symmetry ensures the probability¬ 
preserving evolution of a state, which involves only one 
sub-band or two sub-bands with different k. In this pa¬ 
per we aim at investigating the dynamical behaviors in 
the same model but in the presence of a time-varying 
flux. We have determined that the law of probability 
preservation still holds in the presence of a constant flux. 
When the flux changes adiabatically, we will show rigor¬ 
ously that the Berry phases of all the eigenstates within 
a sub-band are identical complex numbers, which depend 
on the combination of the system parameters, including 
the lattice distortion, imaginary potential, and flux. The 
imaginary Berry phase leads to the amplification and at¬ 
tenuation of the amplitude of an evolved quantum state. 
In contrast to the review [2^ considered in the situation 
where the Hamiltonian has singularities, i.e., diabolic or 
exceptional points, and all the Berry phases in the re¬ 
view are encircling these special points, we would like to 
emphasize the phases investigated in our paper are not 
involving these special points. We also provide some illus¬ 
trative simulations to show that the amplitude control of 
a wave packet can be accomplished by a quasi-adiabatic 
process within a short time. 


This paper is organized as follows. In Sec. m we in¬ 
troduce the non-Hermitian Rice-Mele model and its ex¬ 
act solution. Section Ell is dedicated to show how the 
amplitude modulation is determined by the imaginary 
Berry phase under the time-dependent Hamiltonian. In 
Sec.HVl we investigate wave packet dynamics with its ac¬ 
complishment of amplitude control by a quasi-adiabatic 
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FIG. 1. (Color online) Schematic illustration of the non- 
Hermitian Rice-Mele model driven by a time-dependent ex¬ 
ternal field. It is PT-symmetric with respect to the OO' 
axis in the absence of the field. Constant field breaks the 
PT-symmetry, but keeps the translational symmetry. A time- 
varying field (t) induces eddy field in a direction indicated 
by the red arrow. Together with the distortion, the imagi¬ 
nary potentials can also break the left-right chiral symmetry, 
inducing a direction of the system indicated by the blue ar¬ 
row. In the case that two arrows are either the same or the 
opposite, the dynamics of a state exhibits different behaviors. 


process within a relatively short time. Finally, we give a 
summary and discussion in Sec. El 


II. MODEL AND SOLUTIONS 


We consider a non-Hermitian Rice-Mele model [1^ 
with a flux, which can be described by the following 
Hamiltonian 


2N 

H = + (-1)^' (e*'^c]c,+i + H.c.) (1) 


i=i 


+J{n + iiy)J2i-^y 


on a 2iV-site lattice, where cj is the creation operator of 
a boson (or a fermion) at the jth site with the periodic 
boundary condition C 2 Ar+i = ci. Here the hopping am¬ 
plitude is modulated by the dimerization factor 6, and 
$ = 2N(p is the magnetic flux threading the ring. Stag¬ 
gered complex potential induces the non-Hermiticity of 
the model. 

There are three elements in the structure of the model, 
lattice distortion, imaginary potentials, and flux. As is 
shown in Fig. [TJ the imaginary potentials and distortion 
can break the left-right chiral symmetry (2^ . In addition, 
a time-varying field $ (f) induces eddy field in another di¬ 
rection. The dynamics of a state should exhibit different 
behaviors with different configurations. It turns out that 


imaginar y p otentials can appear in open physical sys¬ 
tems (SOl - l^ . In experiments, the effective magnetic flux 
threading a ring can be realized by rotating the lattice 


In the absence of the imaginary potential and flux, 
this model has been adequately studied in various per¬ 
spectives Q. In the case of purely imaginary potential 
and (j) = 0, the model has PT symmetry, and the dy¬ 
namics has been systematically investigated in the 
frameworks of biorthogonal and Dirac inner products. In 
this paper, we refer particle probability to Dirac proba¬ 
bility. Moreover, we introduce the magnetic flux to the 
Rice-Mele model, which has been employed to control the 
dynamics of wave packet [STl - fs^ . Although non-zero (j) 
breaks the VT symmetry, we will show that this model 
can have full real spectrum. We note that the Hamilto¬ 
nian is invariant via a translational transformation, i.e., 
[T 2 ,H] = 0, where T 2 is the shift operator that defined 
as 


T2^c]T2 = cI+2- (2) 

This allows the invariant subspace spanned by the eigen¬ 
vector of operator T 2 . The single-particle eigenvector of 
T 2 can be expressed as a\ |0) and b\ |0), where 

d = (3) 

= (“) 

^ 3 

satisfying 

T-^alT2 = = e-^%. (5) 

Here, a\. and b\, are two kinds of creation operators of 
bosons (or fermions), with k = 2TrnlN (n S [l,iV]), rep¬ 
resenting the particles in odd and even sublattices. Then 
the original Hamiltonian H can be expressed as 

H = Y.Hk, ( 6 ) 

k 


where 


and 


Hk 


A (fc, (j)) Ja\bk + H.c. 

- (/r -I- iv) J (^a\ak - b\b]}j , 


^ ^) = —E (1 - X5) eWfc/2+0)_ 

A=± 


It is easy to check that 




( 7 ) 

( 8 ) 

( 9 ) 


which ensures us to arrive at the solution in each invari¬ 
ant subspace. 
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Considering the single-particle solution, we 
duce the pseudo-spin operators 

can intro- 

4 = (4)^ = 4 = ^ 

1 , ( 10 ) 

which obey 


[ 4 , 4 ] = 24 , 4 , 4 ] =±4- 

( 11 ) 

Accordingly, has the form 


Hk — ' ^ki 

( 12 ) 


where au is Pauli matrix. Components of the field Bk in 
the rectangular coordinates are 

B^l J = — (1 — 5) cos 4> — (1 + cos {k + (f)), (13) 

Bl/J = (1 - (5) sincj)- (1 + 6) sin (k + (j)) , (14) 

-(/i + w), (15) 

where 

cos6»fc = ^.tanvJfc = (16) 

rij. 

and the field magnitude is 

B, = [iB^,f + {Byf + {B^,fY/^. 

Obviously, 9k can be a complex number even in the case 
with real Bk- 

The eigenvalues of Hk are 

4 = ±Bk (17) 

= ±2J[(^ +w)V4 + ^^ 

+ (l - (5^) cos^ (fc/2 + ()))] 

which give the spectrum of the whole system when all 
possible k are taken. Moreover, the eigenstates of a non- 
Hermitian Hamiltonian can construct a set of biorthogo- 
nal bases in association with the eigenstates of its Hermi- 
tian conjugate. For the present system, eigenstates |'!/'+)’ 

\tjj-) of Hk and |4)’ \V-) are the biorthogonal 

bases of the single-particle invariant subspace, which are 
explicitly expressed as. 



14) = 

| 4 ) = 


— Sin 


Ok 


cos 


— sin ^ 
cos 


(18) 

It is ready to check that the biorthogonal bases 
{\ip\) ,\v\)} (A = ±) obey the biorthogonal and com¬ 


pleteness conditions 


These properties are independent of the reality of the 
spectrum. In this paper, we focus on the case with full 
real spectrum. This happens when (i) z/ = 0, the Hamil¬ 
tonian goes back to a Hermitian system, (ii) /i = 0 and 
v < 26. The phase diagram is true for arbitrary value of 
constant (j), which is the basis for the subsequent investi¬ 
gation of the dynamics with time-dependent flux. 


III. AMPLITUDE MODULATION BY 
IMAGINARY BERRY PHASE 

The aim of this paper is to investigate the effect of the 
time-varying flux on the dynamics of the non-Hermitian 
system. We start with an adiabatic evolution, in which 
an initial eigenstate evolves into the instantaneous eigen¬ 
state of the time-dependent Hamiltonian. 

From Eq. o, we know that iJ is a periodic function 
of (f> with H {(j)) = H {(f> + 2 tt) and Hk {(j)) = Hk (</> -I- 27r). 
Considering the time-dependent flux (f> {t), any eigenstate 
|V'a ( 0 )) return back to \tf}^ ( 0 )) if (j){t) varies adia- 
batically from 0 to 27r, and the evolved state is the instan¬ 
taneous eigenstate More explicitly, we regard 

the flux as a linear function of time, that is, 4> = /3t. And 
the adiabatic evolution of the initial eigenstate \ijjx ( 0 )) 
under the time-dependent Hamiltonian H {4>{t)) can be 
expressed as 

4 ( 0 )) ( 20 ) 

Here the dynamics phase defined by {(j>) and adiabatic 
phase 7 ^ {(j)) have the form 

4(4 = -^^ 44 ) 4 , ( 21 ) 

lk{(t>)=^[ (4 ( 41414 ( 4 ) 4 ( 22 ) 

Jo 

__ 4 2(54d0 

Jo £x [^A - (m + *J^)] ' 

In this paper, we focus on the system with full real 
spectrum. Then the dynamics phase depends only on 
the instantaneous dispersion relation and is always real. 
We have an interest in the adiabatic phase since it may 
have an extra contribution to an evolved state. In the 
following, we will present several features of 7 ^ {(p), based 
on Eq. (1^ . 

It is obvious that we have 74'/’) = ^ 1^6 dimer¬ 

ization factor (5 = 0 , which is a necessary condition for 
nonzero adiabatic phase. In contrast, we have nonzero 
7 ^ {(/)) in the case of zero staggered potential, which in¬ 
dicates the significance of the Peierls distortion to the 
adiabatic phase. It is easy to check that 


4 a ( 4 ) = T exp 


—i H (t) dt 
Jo 


(4 |4) = <5Av4fe', ^14) (41 = 1 - (19) 

A,fc 


44 


7 fc (</>) oc sin {k + 2(1)) 


(23) 
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FIG. 2. (Color online) Plots of the reduced dynamic phase 0^/3 and the geometric phase 7 ^ as functions of 4>, expressed in Eq. 
JUJ and (1221) . for different values of fc in a typical (a) Hermitian system with <5 = —0.15, = 0.05, and v = Q\ (b) non-Hermitian 

system with 5 = —0.15, jj, — 0, and v — 0.05. Both systems have full real spectrum. The real (imaginary) parts of the phases 
are plotted in red (black), (a) It is shown that both phases are real for the Hermitian system, while the adiabatic phase is 
complex for the non-Hermitian system. The dynamic phases and the real part of adiabatic phases in both cases have a slight 
difference. The shapes of all the phases are in agreement with the approximate expressions in Eqs. (O, dSl, ISzl), and ([381). 
Here the reduced dynamic phase 0^/3 and the geometric phase 7 ^ are expressed in units of radian. 


which leads to 9^7^ {4>) = 0 at (j>c = it 12 — k 12. This 
fact indicates that both the real and imaginary parts 
of 7 ^ {(j>) experience a maximal (or minimal) change at 
this point. Furthermore, including Hermitian and non- 
Hermitian systems, we note that 

dlt jmr) da^ ( titt ) _ ^ 

dk dk ' ^ ’ 

In view of the fact that for an arbitrary function 
g (cos {k -h 2(j))), one can always obtain 

g (cos (fc -I- 2(j))) dcj) 

= ^5(cos(fc-h2(/)))”’'= 0. (25) 

This shows that the dynamic and Berry phases are k- 
independent after (f> varying mr, which ensures that any 
arbitrary initial state involved in the upper or lower band 
solely can revive back exactly after (j) varying 27r. 

For a Hermitian system, the adiabatic phase is always 
real, which ensures the probability preserving evolution. 



While, the probability of an evolved state changes due 
to the imaginary part of the adiabatic phase in a non- 
Hermitian system. The gain or loss of probability de¬ 
pends on the sign of the imaginary phase. There are 
several rigorous results for the phase. 

For given parameters the adiabatic phase of 

an eigenstate in A band obeys the identity 


Ik {mr) = -7fc {-mr) , (26) 

which is obtained from Eq. (1221) . owing to the fact 

£7 (</> - rnr) = el ((/»). (27) 

This means that the direction of </> can determine the sign 
of Imy^ {mr), controlling the probability of an evolved 
state. 

In addition, the sign of Rey^ {(j)) and Imy^ {(j)) could 
also depend on the sign of 6 and v via the Eqs. (1251) and 
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(l29l) . that is, 


Re 7 fe (</>) = sgn (5) 


-2 |(5| J‘^d(t> 


{Bkf + 


(28) 


Im 7 fc ((?i) = sgn((5i/A) 


-2J^\5v\d(l) 


Bk 


{Buf + J^iy^ 


which yield 


bk (<(>)] 5 = - [7fc 

bfe W], = [7fe (</>)]!,• 


(29) 


(30) 

(31) 


Together with the Eq. (EH), we show that the sign of 
Im ( 7 ^) is determined by the following expression that 

sgn [imy^ (A'nTr)] = —sgn (^i5AA') (A'= ±), (32) 

which directly results in the amplification or attenuation 
of an evolved state. 

Besides the exact results, it is useful to arrive at the 
whole profile of phases as a function of a group of pa¬ 
rameters S, fi,v, X}. The approximate expressions 
of the phases for Hermitian and non-Hermitian systems 
could be obtained from straightforward derivations as fol¬ 
lowing, respectively. 

For a Hermitian system (j^ = 0), we have 


7fc {<!>) 


sgn ((5/2) 


{tan ^ Qk (0) — tan ^ 0^ (^) (33) 


-sgn(^A) [tan"^ (ImITa: (/>)) - tan“^ (iMlTfe (0))]}, 
and 

A 


(</>) 


-{Bk (0) (fc - tt) - Bk (</>) {k-\-2(j)- tt) 


J (4^^ + Bk (()>) - 2 1^1 JQk (</) 

Bk ( 0 ) - 2 | 5 | JQk ( 0 ) 


}, 


(34) 


where 


Qk (</>) = \/(5 2 - 1 {k + 2(j)- tt) /2, 

T/c (</>) = J0fc (</) /Bk {(/) , 


(35) 

(36) 

are even functions of (5, /x, and A. It is shown that 7 ^ {(/) 
is a real number, preserving the probability. 

While, for the non-Hermitian system (/x = 0, (/ < 25), 
we have 


Ik 


sgn ((5/2) 


(tan ^ Qk (0) — tan ^ Qk {(j)) (37) 


-isgn (i/A) [tanh“^ (|ix| Tk ((/)) - tanh“^ (|i/| Tfe (0))]}, 
and 

A 


W - 4 ^ 

J (4^2 _ ^ 2 ^ _ 2 1^1 J 0 ^ 


(Bfc (0) {k - tt) - Bk {(/) {k + 2(l)- tt) 

}, 


Bk ( 0 ) - 2 |(5| JQk ( 0 ) 



(a) 



(b) ^{t)ln 

FIG. 3. (Color online) Numerical simulations for the time 
evolutions with different /?. Amplification factor A/. {(/) and 
fidelity // {(f)) for the initial state with k = 7r/25 and A = 
— are plotted for the system with S = 0.15, fj. = 0, d = 
—0.2, and N = 50. We can see the evolution results are 
close to that, A~^^^ {tt) = 27.58, f ~^25 (’^) ~ obtained in 
the adiabatic limit. This indicates that the amplitude control 
could be realized with a high fidelity via a quasi-adiabatic 
process. 


where 0 ^ {(/>) and Tk {(f) have the same form as above, 
but even functions of 5, v, and A. The remarkable feature 
of 7 ^ is that it is a complex number. The approximate 
expression shows that both real and imaginary parts of 
7 ^ {(j)) are flat functions of (p except for the region around 
the point (pc-, in which they experience a drastic change. 
The key feature of an imaginary phase is its sign, which 
affects the amplitude of the evolved eigenstate directly, 
determining the gain or loss of the probability. 

To verify and demonstrate the above analysis, numeri¬ 
cal simulations are performed to investigate the dynamics 
behavior of a quasi-adiabatic process. We compute the 
time evolution of an eigenstate by using a uniform mesh 
in the time discretization for the time-dependent Hamil¬ 
tonian H{t). The amplification factor (gain) is defined 
as 


A^k W = ||V'!(0))| 


T exp 


H {t) dt 


bx (o)> 


(39) 


(38) 
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FIG. 4. (Color online) The comparison between the disper¬ 
sion relation (blue dot) and the trajectory of the center of the 
wave packet (color contour map). The simulation is computed 
as follows: the wave packet with fco = 7r/2, a = 0.05, and 
Na = 1900, in the lower band of the system with 5 = —0.15, 
^ = 0, u = 0.05, and N = 1000, subjected in the external 
field with (3 = 0.001. This shows that the two are in close 
proximity owing to the fact that the contribution of the geo¬ 
metric phase to the wave packet position becomes negligible 
when compared with the dynamic phase. Only the imaginary 
part of the geometric phase takes an important role in the 
dynamics of the wave packet. 


which is the ratio of the output magnitude to the input 
magnitude of an eigenstate. We use the fidelity {(/>), 
which is defined as 


/fc (<(>) = 

(^A {.(t>)\T exp 

—ifn {t) dt 

I^A ( 0 )) 



Jo 



to describe the derivation between adiabatic and quasi- 
adiabatic processes. For an adiabatic process (,5 —>■ 0), 
we have A^{4>) = exp [-Im ( 7 ^ (i^))] and = 1 . 

We can employ (tt) to describe the amplification fac¬ 
tor for an arbitrary quantum state. The computation 
is performed by using a uniform mesh in the time dis¬ 
cretization for the time-dependent Hamiltonian H{t). As 
an example, in Fig. [3l we show the evolution of A^ ((/>) 
and {(j)) for different values of (3. The plot in (a) and 
(b) shows the quasi-adiabatic process can be close to the 
adiabatic one. 


IV. WAVE PACKET DYNAMICS 

Before starting the investigation of the wave packet 
dynamics in the present system, we would like to give 
a brief review on the dynamics in a uniform ring. In 


the case of 5 = ^ = iz = 0 , it has been shown that the 
dynamics of a wave packet is the same as that driven by 
a linear field with strength /3, according to the quantum 
Faraday’s law [s^. Furthermore, it is turned out that 
the center path of a wave packet driven by a linear field 
accords with the dispersion of the Hamiltonian in the 
absence of the field within the adiabatic regime [i^ . 

Xc {(j)) = Xc (0) -b i {e'='= {(j)) - (0)} , (41) 

where {(j>) is the dispersion relation and kc is the cen¬ 
tral momentum of the wave packet. 

Now, we switch gears to the case of the present model. 
We note that the geometrical phase takes a role in the 
dynamics and the dynamics of the wave packet cannot 
be simply understood in terms of a semiclassical picture 
[dll - li^ . Notice that the trajectory of a wave packet is 
essentially not only determined by the dispersion relation 
for the field-free system but also by the geometric phase. 

In other words, the dependence of the geometric phase 
on k should be taken into account. 

However, in the adiabatic limit, /3 is very small. The 
contribution of the geometric phase to the trajectory be¬ 
comes negligible, that is, 

{4^) « al (4 ). (42) 

The main effect of the geometric phase on the dynamics 
of the wave packet is then just the modulation of the 
amplitude. 

To verify and demonstrate the above analysis, numer¬ 
ical simulations are performed to investigate the dynam¬ 
ics behavior. We compute the time evolution of the wave 
packet by the same method as mentioned above. The 
initial Gaussian wave packet has the form 

. . 2N 

G^-(O)) = \j) (43) 

with the velocity ko G [0, 27r], centered at the VAth site, 
where |j) = cj | 0 ) and Hi is a normalization factor. 

For one thing, we consider the trajectories of the wave 
packet with different /3, as comparison to the dispersion 
relation. From the plots in Fig. 01 we find that for small 
/3, the trajectory accords with the dispersion well, while 
as (3 increases, the deviation becomes obvious. 

For another thing, we investigate the flux-controlled 
probability of the wave packet. It can be rewritten in the 
form 

Gl^{0))=Y.{9X\4X)+9t\4t)). (44) 

k 

Here, we do not give the explicit expression of the coef¬ 
ficient since the following analysis is independent of 
g^. Through an adiabatic evolution, we have 

G^"(2W/9)) = E \i,X) + |^^)) , 

k 

( 45 ) 
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FIG. 5. (Color online) The profiles of the time evolution of a wave packet in several typical cases. The initial wave packet is in 
the form of Eq. (143 II with fco = 1.47r, a = 0.05, and Na = 1000, and the external field varies with = 0.001. (a) In the case of 
S = —0.7, n — 0, and v = 1.3, the plot is shown that the probability of the sub-wave packet in one band increases but the one 
in another band decreases. In the cases of (b) and (c), we take the same parameters with (a) bnt opposite (p and S, respectively. 
Comparing with the profile in (a), we can see the completely opposite behaviors in (b) and (c). This is in agreement with our 
analytical prediction. 






(a) 


f-T 


(b) 


t-T 


(c) 


t-x 


FIG. 6. (Color online) Profiles of the reduced energy defined in Eq. (1481) of an evolved wave packet in Eq. (1431) is controlled by 
a Ganssian-type time-varying field in Eq. (|43). The parameters of the wave packet and system are ko — l.bvr, <5 = —0.7, /r = 0, 
and ly = 1.3. (a) The shapes of the field as a fnnction of time in Eq. ()47|l with different a. (b) Plots of E{t) of the field in (a), 
(c) Same as (b) but with the opposite flux. The simulation resnlt is in agreement with the result E{oo) = —1.14 for (b) and 

I. 14 for (c), obtained in the adiabatical limit. This indicates that the amplitude of the wave packet can be well controlled via 
a quasi-adiabatic process in a short time. Here the time t — r and the reduced energy E (t) are expressed in units of 1/J and 

J, respectively. 


where is a real number, and = —Im 7 ^ (27rn). In 
the case of ^ 1, we have 

Gl-{2nn/P)) ^ IV'a) • (46) 

k 

It is shown that only the component in either upper or 
lower sub-band survives. It also seems that the final state 
collapses to one of two sub-bands in the context of Dirac 
probability. The sign of C is crucial for the direction of the 
collapse. We compute the evolution for two cases with 
opposite flux. The result plotted in Fig. [S]shows that the 
evolved wave packet in upper or lower sub-band survives 
for two different varying fluxes, which is in agreement 
with our prediction. 

In practice, the flux control can be implemented by 
a pulsed flux. We simulate this process by a Gaussian- 


shaped flux with the form 



which contributes 27r flux during the process. To charac¬ 
terize the feature, we introduce the reduced energy 


E{t) 


J2k,\ 





(48) 


And in the adiabatic limit, we have A (oo) ~ or , 
i.e., it converges to a positive (negative) constant if the 
upper (lower) sub-wave packet survives. Alternatively, 
we can replace by A in Eq. (H51) to redefine E (t), which 
leads to E (oo) Ri 4-1 or —1 even for a non-adiabatic 

































process. Here we take the former, because E (oo) can 
indicate the deviation between the adiabatic and non- 
adiabatic processes. In Fig. El results are plotted for 
different values of ct, which determine the rate of the flux 
change. These results clearly demonstrate the finding of 
this paper that the external field can be utilized to con¬ 
trol the quantum state on demand via a quasi-adiabatic 
process within a relatively short time. 

V. SUMMARY 

In this paper, the particle dynamics of the non- 
Hermitian Rice-Mele model driven by a time-dependent 
external field has been theoretically investigated. The 
analysis shows that the Berry phase can be a complex 
number in the non-Hermitian regime. This results in 
the amplification and attenuation of the amplitude of an 
evolved quantum state. We found that it can have full 
real spectrum for any constant field, and the Berry phase 
with respect to a varying field has a constant imaginary 


part for an arbitrary initial state either in the upper or 
lower energy sub-band. The dependence of the imagi¬ 
nary part of the Berry phase on the parameters, such as 
lattice distortion, imaginary potential, and the direction 
of the flux, was explicitly presented. Numerical simula¬ 
tion indicates that the amplitude control of a wave packet 
can be accomplished by a quasi-adiabatic process within 
a relatively short time. 


ACKNOWLEDGMENTS 

We acknowledge the support of the National Basic Re¬ 
search Program (973 Program) of China under Grant No. 
2012CB921900 and CNSF (Grant No. 11374163). 


REFERENCES 


[1] S. Pancharatnam, Proc. Indian Acad. Sci. A 44 247 
(1956). 

[2] M. V. Berry, Proc. R. Soc. London A 392, 45 (1984). 

[3] D. Xiao, M. -C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 
1959 (2010). 

[4] J. C. Garrison and E. M. Wright, Phys. Lett. A 128, 177 
(1988). 

[5] G. Dattoli, R.Mignani, and A. Torre, J. Phys. A 23, 5795 
(1990). 

[6] C. Z. Ning and H. Haken, Phys. Rev. Lett. 68, 2109 
(1992). 

[7] D. J. Moore and G. E. Stedman, Phys. Rev. A 45, 513 
(1992). 

[8] Xiao-Chun Gao, Jing-Bo Xu, and Tie-Zheng Qian, Phys. 
Rev. A 46, 3626 (1992). 

[9] Marcel Pont, R. M. Potvliege, Robin Shakeshaft, and 
Philip H. G. Smith, Phys. Rev. A 46, 555 (1992). 

[10] S. Massar, Phys. Rev. A 54, 4770 (1996). 

[11] Y. C. Ge and M. S. Child, Phys. Rev. A 58, 872 (1998). 

[12] Robert S. Whitney, Yuriy Makhlin, Alexander Shnirman, 
and Yuval Gefen, Phys. Rev. Lett. 94, 070407 (2005). 

[13] H. Mehri-Dehnavi and A. Mostafazadeh, J. Math. Phys. 
49, 082105 (2008). 

[14] Alexander 1. Nesterov and S. G. Ovchinnikov, Phys. Rev. 
E 78, 015202(R) (2008). 

[15] Xiao-Dong Cui and Yujun Zheng, Phys. Rev. A 86, 
064104 (2012). 

[16] Shi-Dong Liang and Guang-Yao Huang, Phys. Rev. A 
87, 012118 (2013). 

[17] C. Li, G. Zhang, X. Z. Zhang, and Z. Song, Phys. Rev. 
A 90, 012103 (2014). 

[18] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 
5243 (1998). 

[19] C. M. Bender, D. C. Brody, H. F. Jones, and B. K. Meis- 
ter, Phys. Rev. Lett. 98, 040403 (2007). 


[20] A. Mostafazadeh, Phys. Rev. Lett. 102, 220402 (2009); 
Phys. Rev. A 80, 032711 (2009). 

[21] S. Longhi, Phys. Rev. B 80, 165125 (2009). 

[22] X. Z. Zhang, L. Jin, and Z. Song, Phys. Rev. A 87, 
042118 (2013). 

[23] X. Q. Li, X. Z. Zhang, G. Zhang, and Z. Song, 
arXiv: 1409.0420 

[24] G. Zhang, X. Q. Li, X. Z. Zhang, and Z. Song, 
arXiv: 1412.7630 

[25] T. E. Lee, F. Reiter, and N. Moiseyev, Phys. Rev. Lett. 
113, 250401 (2014). 

[26] S. Longhi, Phys. Rev. A 88, 062112 (2013) 

[27] W. H. Hu, L. Jin, Y. Li, and Z. Song, Phys. Rev. A 86, 
042110 (2012). 

[28] Ingrid Rotter, J. Phys. A: Math. Theor. 42 153001 
(2009). 

[29] M. J. Rice and E. J. Mele, Phys. Rev. Lett. 49, 1455 
(1982). 

[30] J. Dalibard, Y. Castin, and K. Mplmer, Phys. Rev. Lett. 
68, 580 (1992). 

[31] K. Mplmer, Y. Castin, and J. Dalibard, J. Opt. Soc. Am. 
B 10, 524 (1993). 

[32] R. Dum, P. Zoller, and H. Ritsch, Phys. Rev. A 45, 4879 
(1992). 

[33] M. B. Plenio and P. L. Knight, Rev. Mod. Phys. 70, 101 
(1998). 

[34] A. J. Daley, Adv. Phys. 63, 77 (2014). 

[35] S. Tung, V. Schweikhard, and E. A. Cornell, Phys. Rev. 
Lett. 97, 240402 (2006). 

[36] R. A. Williams, S. Al-Assam, and C. J. Foot, Phys. Rev. 
Lett. 104, 050404 (2010). 

[37] T. J. Osborne and N. Linden, Phys. Rev. A 69, 052315 
(2004). 

[38] S. Yang, Z. Song, and C.P. Sun, Phys. Rev. A 73, 022317 
(2006). 



9 


[39] W. H. Hu, L. Jin, and Z. Song, Quantum Inf. Process 
12(11), 3569 (2013). 

[40] S. Lin, X. Z. Zhang, and Z. Song, Phys. Rev. A 90, 
063411 (2014). 


[41] F. Bloch, Z. Phys. 52, 555 (1929). 

[42] N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Harcourt, Fort Worth, TX, 1976), Chap. 12. 

[43] C. Kittel, Quantum Theory of Solids, 2nd ed. (Wiley, 
New York, 1987). 



